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Abstract 

The Lie algebraic integrability test is applied to the problem of classification of inte- 
grable Klein-Gordon type equations on quad-graphs. The list of equations passing the test 
is presented containing several well-known integrable models. A new integrable example 
is found, its higher symmetry is presented. 
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1 Introduction 

We study the integrability problem for the quad-graph equation of the form 

= (1) 

Such kind equations have a large variety of applications in physics, biology, architecture, etc. 
Here the unknown u = u(m,n) is a function of two discrete variables m, n. For the sake 
of convenience we use the following notations: Uk = u(m + k, n), Uk = u(m,n + k), u^i = 
u(m + 1, n + 1). Function / is supposed to be locally analytic, it depends essentially on all 
three arguments. In other words equation (TTj can be rewritten in any of the following forms 

u i,j = f^i^UuUj), with i = ±1, j = ±1. (2) 

Nowadays various approaches are known for studying integrable discrete phenomena. The 
property of consistency around a cube [TJ, has been proposed as the integrability criterion 
for quadrilateral difference equations [2], [3J, [I]. Symmetry approach to the classification of 
integrable systems is adopted to discrete case [5J, [B], [TJ, [S], [H]- Another characteristic property 
of an integrable equation is the vanishing of its algebraic entropy [TU] . Alternative methods are 
used in [11], [12], [13]. In this article we study discrete phenomena from some different point 
of view. 

Years ago it was observed that characteristic Lie algebras, introduced in [TJ] , in the case of 
integrable hyperbolic type PDE's like sine-Gordon and Tzitzeica-Zhiber-Shabat equations have 
a very specific property. The dimensions of the linear spaces spanned by multiple commutators 
of the generators grow essentially slower than in generic case. In [15] the problem of rigorous 
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formalization of the notion of "slow growth" has been discussed. A conjecture was suggested 
and checked by applying to classify integrable equations of the form u XjV = f(u,u x ). 

In the recent article [16] we introduced and successfully tested a classification scheme (called 
algebraic test) based on investigation of multiple commutators of characteristic vector fields 
defined by equation ([T]). Now we use the algebraic test (an explanation is given in Section 2 
below) to classify discrete Klein-Gordon type equation on quad graph 

ui,i +u = g{u x (3) 

The list of equations which partially passed the algebraic test is studied additionally by applying 
the symmetry test. It is remarkable that the final list contains in addition to the well-known 
equations as the discrete potential KdV equation, the discrete Liuoville equation and a discrete 
analogue of the sine-Gordon equation also a new integrable model (see below the discussion on 
Corollary of theorem 2 and theorem 3): 

ui,i = — ■ 4 

U\OL + U\Ui) 

Remind that the problem of complete description of integrable Klein -Gordon equations 

Ux,y = g(u) (5) 

has been solved years ago by A.V.Zhiber and A.B.Shabat [17]. The authors proved that the 
only integrable equations in (jSJ) are 

linear equation; 

the Liouville equation u XyV = e u ; 

the sine-Gordon equation u X}V = siiau; 

the Tzitzeica-Zhiber-Shabat equation u x>y = e u + e~ 2u . 

The article is organized as follows. In section 2 we define characteristic vector-fields and the 
test Lie ring, formulate the algebraic test, consider an illustrative example. The classification 
problem for the model of the form ui } i + u = g{u\ + U\) passing the test is investigated in 
section 3. The result of classification is summarized in Theorems 3 and 5. 



2 Algebraic test and classification scheme 

Note that for any integers p and q the variable u Ihq is expressed in terms of the variables 
{v,i, Uj}°°j = _ 00 in a recurrent way. Hence the variables Uj, Uj are called dynamical variables. 
They considered as independent ones. For example, 

v>2,i = f(u 1 ,u 2 ,u 1)1 ) = f(u 1 ,u 2 ,f{u,u 1 ,u 1 )), 

and 

U-l-2 = / _1 '~ 1 ('U-l,W_l-l,'U-2) = /~ 1 '" 1 (W-l,/" 1 '~ 1 (w,'ti-l,«-i),«-2), 

Introduce the shift operators D and D, shifting the first and respectively the second 
integer arguments: Dh(m,n) = him, + l,n), Dh(m,n) = h(m,n + 1). Explain the ac- 
tion of the operators D and D on the functions of dynamical variables. For the function 
h = h(ui,Ui_i...,Ui',Uj,Uj + x, ...,%') we have 

hk = D k h = h(iii + k, Ui + k-i---,Ui' + k, Ukj, Ukj+i, ...,Uk,f) 
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and similarly 

hk = D k h = h(ui : k, Ui > k-i---,Ui> 7 k, Uk+j, Uk+j+i, ...,Uk+j>). 

In these two formulas one has to replace the double shifts u a> p through dynamical variables as 
it is discussed above. 

Characteristic vector fields for the hyperbolic type PDE were introduced by E.Goursat 
[18] in 1899. They provide a very effective tool for classification of Liouville type integrable 
systems. For instance in [18] a list (almost complete) of hyperbolic type PDE was found ad- 
mitting integrals in both directions by using a method based on the notion of characteristic 
vector fields. Interest to the subject renewed after the paper [TJ], where the exhaustive de- 
scription of exponential type systems is given admitting complete set of integrals. The concept 
of characteristic vector fields was adopted to the quad graph equations in [T9]. In the recent 
papers [20], [2TJ, [22J the characteristic vector fields were used for the purpose of classification of 
semi-discrete Liouville type equations. All of these studies concern with the Liouville type inte- 
grability. In [15] and [16] some new applications of these important notions are suggested. Let 
us remind some necessary definitions. At first suppose that equation ([1]) admits a non-trivial 
n-integral, i.e. there is a function / depending on a finite number of the dynamical variables 
/ = I(u^j,u^j + i, ...Uk) satisfying the equation DI = I. In an enlarged form the last equation 
is 

I( r -j+h r -j+2, r, u, f, fi, ...fk-i) = I{u-j, M-j+i, ...u k ), (6) 

where the function r = u~i, u\) is defined in ([2]). Since the right hand side of (jSJ) does 

not depend on the variable U\ we find -g^DI = -rj^I = 0. It implies the equation YI = where 
the operator Y is defined as follows Y := D~ l -^-D. Direct computations give (see [16], [T9] ) 

Tr d d 1 d did 

Y = — + x— + + xx\— h + ..., 7 

where x = ^-i ^/K^ih = _ d^-\u,^,u- 1 )/du We u th operators x := J- and Y 

characteristic vector fields. Now it is evident that the map f(u,Ui,Ui) =^ Y is correctly 
defined for any / due to the formula (jTJ). 

Denote through T the set of vector fields obtained by taking all possible multiple commuta- 
tors of the operators X and Y, and taking all linear combinations with coefficients depending 
on a finite number of the dynamical variables U-i, u, u±±, u±2, ■■■ ■ Evidently the set T has 
a structure of the Lie ring. We call it test ring of the equation ([1]) in the direction of n. In a 
similar way one can define the test ring T in the direction of m. 

Notice that for Liouville type integrable equations of the form (JTJ) both rings T and T are 
of finite dimension. Actually the test ring is a subset of the characteristic Lie ring [19j-|22j. 

Denote through Vj the linear space over the field of locally analytic functions spanned by 
X, Y and all multiple commutators of X and Y of order less than or equal to j such that: 

V = {X,Y}, V 1 = {X,Y,[X,Y}}, .... 

Remind that such kind sequences of linear spaces have important applications in geometr}{§ 
(see [23], [21]). Introduce the function A(k) = dim\4 + i — dimVfc. The following conjecture is 
approved by numerous examples. 

Conjecture (algebraic test). Any integrable model of the form (Q]) satisfies the following 
condition: there is a sequence of natural numbers {tk}^^ such that A(tk) < 1. 
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Ring T admits an automorphism, generated by the shift operator D, 

T3 DZD- 1 e T, (8) 

which plays crucial role in our further considerations. It is important that X and Y considered 
as operators on the set of functions depending on the variables m_i, u, u±i, u± 2 , ... satisfy the 
following conjugation relations [16] 

DXD- 1 = pX and DYD~ X = —Y, (9) 

x 

where p = D(— J^'" -1 '" -1 ^ ) = df\-i(uulu {Um r m deed, specify the coefficients of the 

operator DXD' 1 = J2 a i^r + Po^~[ by applying it to the dynamical variables and find that 
a,i = DXD~ 1 u i = for any integer i. Moreover p = DXD~ l u_\ = Dl/ _lrl (u, w_i, m_i) = 
D( a/ " 1 '" 1 ^'"- 1 ' fi - l) ). In a similar way one can prove the second formula. Really apply the 
operator DYD~ l = + d^—^ to ui and find Cj = D(Yuj-\) = Yuj. Then evaluate d = 

DY D~ x u~\ = fu 1 ' 1 + —fulf 1 - Since m_i_i = m_i, m_i) one gets the equation u = 

Mi, Mi), Mi, Mi). Let us differentiate it with respect to u\ and find DD( 9 ^-^ — )m[ + 
DD( d f~^' ' * ) = or the same + sf 9 ^ du ' = 0. Now due to the equation 

x = Z) _1 (J^) one concludes that d = 0. 

Lemma 1. Suppose that Z = Y^oo^j^T T satisfies the following two conditions: 

1) DZD^ 1 = cZ for some function c; 

2) bj = for some fixed value of j = j . 

Then Z = 0. 

Proof of the Lemma can be found in [16J. 

Example. As an illustrative example we consider the following discrete Liouville type 
equation (see [25] ) 

1*1,1 = -(«! - 1)(W! - 1)- (10) 
M 

Find first an explicit form of the characteristic vector field Y. Since / = -{v,\ — 1)(mi — 1) and 
§- = - 1) then 



x = j D- 1 -(mi-1) = — (11) 

M M_i 

Express «i,_i through the dynamical variables due to equation (TlUj) . Apply the operator D^ 1 
to both sides of f fTUj) and find Mi = ^— — l)(u— 1). Comparison of the last equation with 
(ITT]) yields x = -r^-, therefore 

<9m m — 1 C?Mi M (9u_i [u — 1)(mi — 1 ' 

For this equation p = Jtl = I, hence DXD' 1 = pX, DYD~ X = pY . Evaluate D[Jf, Fp^ 1 = 
p 2 [X, y] +pX(p)F-py(p)X, :r smce X(p) = Y(p) = we get DLY, = p 2 [X, Y], moreover 

[X, y] = Y^jLi a j~(jt~ Jt ~ a -jdu~' Now due to the Lemma 1 one obtains [X, Y] = 0. So dimension 
of the ring T for equation ( ITU]) equals two and A (A;) = for all k > therefore equation ( ITU]) 
passes the test. 
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3 Equations of the form + u = g(u\ + U\) 

In this section we apply the Conjecture to the following Klein-Gordon type particular class 
of discrete model (0Q) 

ui t i + u = g(u! + Mi). (13) 

where the function g is to be determined. 

Classification scheme. Obviously equation (CQ) passing the algebraic test above should 
satisfy one of the conditions: 

i) A(0) < A max .(0) = 1; 

ii) A(0) = A max {0), A(l) < A max (l) = 2; 

hi) A(0) = A max (0), A(l) = A max (l), A(2) < A max (2) = 3; 

iv) A(0) = A max (0), A(l) = A max (l), A(2) = A max (2) and A(k) < 1 for some k > 2. 

Where A maa; (fc) stands for the greatest value of A(k) for equation (0Q) when f(u,ui,ui) ranges 
the class of arbitrary functions. 

Remark that in the case of equation ({TBI investigation of the first three particular cases 
allows one to extract a very short list of equations expected to be integrable. The list is 
exhaustive because the case iv) is never realized (see below Corollaries of Theorems 2 and 4). 

Introduce vector fields i?i = [X,Y], Pi = [X,^], Q 1 = [Y, Ri], R 2 = [X,Q X ], W = [Y,Q 1 ], 
Z = [X, Pi]. Using these vector fields we can span in addition to Vq and V\ two more linear 
spaces: 

V 2 = V 1 + {P 1 ,Q 1 }, V 3 = V 2 + {W,Z,R 2 }. 

In order to evaluate A(k) we will use the automorphism (jHJ). At first evaluate the factors x 
and p in formula ([9]) for the case (fT3jl . We have x = p = g'(g~ l (ui + u_i)), where the function 
P = g~ 1 (a) is the inverse to the function a = g(/3). Inversely, knowing x = x(u\ + one 
can recover g(/3) by using the equation 

» = = /(^M)'* = / = J £y (") 

Specify the action of the characteristic operators X and Y on the variable x. Evidently, 
Xx = x', Yx = xx 1 . It is found by direct calculation that 

x' 

DRxD- 1 = R x Y- x'X, 

x 

'2 

DP\D~ l = xP\ — x'R\ — rY — xrX, r = x , 

x 

1 ' " 

DQiD' 1 = -Q 1 + -R 1 -—Y-x"X, 

XXX 

DWD- 1 = ±-W + (—-^ T )R 1 -—Y-x"'X, (15) 



/ y» 2 rift / y» 2 



'3 
X ^ 



DZD = x Z + (x — 2xx )R\ — qY — xqX, q = xx — 2x x + 

x 

I '2 / " 

DR 2 D- 1 = R 2 --Q 1 + x'P 1 -—R 1 -sY-xsX, s = x" - — , 

X X X 
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Let us study the set G of all multiple commutators of X and Y. 

Lemma 2. The coefficients of any operator in G are functions of a finite number of the 
dynamical variables x, x±i, x±2, ... ■ 

Proof. Here x(a) is a fixed function of the argument a = u\ + Therefore one can 
write x' = <fr(x) for some function 0. Then X(x) = —<f>(x) and Y{x) = x(f)(x) =: ip(x). By 
using the conjugation relations DXD~ l = xX and DYD^ 1 = -Y one derives that X(xj) = 
4>i(x, xi, ■■■Xj) and Y(xj) = ip^x, x\, ...Xj). Similarly X(x-j) = <p~^(x, X-\, ...X-j) and Y(x-j) = 
i/)~i(x, X-i, ...X-j). Now evidently R± = X(x)-£^ + X{-^-)-^—^ + ^( XXl ^^ + • ■ ■ satisfies the 
statement of the lemma. Due to the formulas Ri(x) = X(x)x' and DR\D~ l = R± + —Y — x'X 
one gets R±( ). Obviously the proof can be completed by using induction. 

Theorem 2. Suppose that an equation of the form U~3\) satisfies one of the conditions 
1) — Hi) of the Classification scheme then function x = x(a) solves the following ordinary 
differential equation 

x' 2 = (x 2 + 1)7 + xv (16) 

with the constant coefficients 7, v. 

Proof of the theorem 2. Begin with the case i), suppose that A(0) = 0, then we have 
Ri = XX + iiY. It is evident that R x = X(x)^ + X = ^ and Y = £ + x-^- + ... hence 
A = fx = and therefore R\ — 0. By applying the automorphism above to both sides of the 
last equation one finds 

-—Y-x'X = 0. 

x 

Since X and Y are linearly independent we get equation x' = which is a particular case of 
f Tl6l) . Evidently its solution is x = c and due to (fT4l) it can be found that (3 = g~ l (a) = -a-\-c\. 
Thus our equation a = g((3) (see ( ITS]) ) is linear U\,\ + u = c{u\ + u\ + c\). For this case 
dimT = 2 so that A(k) = for k > 0. In a similar way one checks that condition ii) leads to 
ffTB"]) . Indeed suppose that A(0) = 1 and A(l) < 2 then we have 

Px = vQ x + eR x . (17) 

Here due to Lemma 2 functions v and e might depend only on x, x±i, x±2 .... Apply the 
automorphism (jSJ) to both sides of (|17p then simplify due to equations (|15|) : 

I ' " ' 
xtuQt + eRi) - x'Rx -rY - xrX = D(y)(-Qx + - — F - x"X) + D(e)(R l Y - x'X). 

rr> ry rr* rr* 

JU JU i\j iAj 

Comparison of the coefficients before linearly independent operators gives rise to the conditions 

Qi : xv = -D(v); 

x 

Rx : -x' + xe = —D(v) + D(e); 

x 

Y : -r = ——Die) + —D(v); 

x x 

X : -xr = -x'D(e) - x"D(v). 

Simple analysis of these equations implies that equation fTlTj) holds if and only if the following 
three conditions valid: v = 0, e = const, x' = e(x — 1). Actually, under these conditions 
the last two equations above are satisfied automatically. In a similar way one can check that 
Qi = vPi — eRi is equivalent to the same three conditions. Hence if A(l) < 2 then A(l) = 
and consequently A(k) = for any natural k > 1. Thus in this case dimT = 3. 
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Let us suppose now that A(0) = 1, A(l) = 2 and A(2) < 2 which corresponds to the case 
Hi). At first consider the case when Z is linearly expressed through the other vector fields in 
subspace V3: 

Z = 7^! + SP X + eQi + <\>R 2 + ij>W. (18) 



Apply the automorphism (jSj) to both sides of equation ffl8|) and compare coefficients before 
linearly independent operators 

W : x 2 ip =D(ip) — 



x 2 



R 2 : x 2 <P =D((j)), 

Qx: x 2 e =D{e)- + -<j), 
x x 

P 1 : x 2 5 = D{S)x + x'cj), 

Ri: x 2 7 +x' 2 -2xx" = D(6)x' + D(j). 

Since x = x{u\ + u~\) we have ip = 0, (ft = 0, e = 0, 5 = 0, 7 = const. Comparison of 
the coefficients of X and Y gives one more equation xq = 72'. Finally we get two ordinary 
differential equations for x which are absolutely the same as in the case of the equation Mi,].— u — 
g(ui - Mi) (see (T6J): 

— 2xx x + x = 7X , (x — 1)7 + x — 2xx = 0. 

The compatibility condition of these equations is equivalent to equation ([TBI) . In this case we 
have Z = jR±. It is remarkable that x solves equation (fToT) if and only if W is linearly expressed 
through the other elements of V3 and then W = jRi- And the last possibility is when R 2 is 
linearly expressed through X, Y, Ri, Pi, Qi, W, Z. In this case x solves the equation x' = 0. 
The proof of the theorem is completed. 

In order to find x = x(a) evaluate the integral 

r dx 

H(x) := I = a - oo. (19) 

(X 2 + 1)7 + XV 



For the case 7 7^ the answer is given by the formula 



H(x) = ln(2 v / x 2 + l + xfe + 2x + 6), b=-. (20) 

>Jl 7 

Now find x by solving the equation H(x) = a — a^. 

x(a) = - e ^ {a - ao) - — _(!_ Zl) e -v^(«-«o)_ 
v ; 4 2 7 v 4 7 2 ; 

In order to get the corresponding quad-graph equation (fi~3l) integrate again: 



Integration gives 



= ^-]n 



e ^t{a-a ) _ ^ + 2 

7 



+ /So- (22) 



Then the equation searched is given by the formula (take y/^y = 1) 

a ie Ul < 1+Ul+ai+u + a 2 e Ul ' 1+u + a 3 e Ul+Ul + a 4 = 0, 

where a\ = e~ ao_/30 , a 2 = — e~ a °, a 3 = (2 — 6)e~^°, a 4 = 2 + b. By setting u = \nv one reduces 
it to a quadrilinear form 

aivi^viviv + a<iV\^v + a%V\V\ + a 4 = 0, (23) 
with arbitrary constants a,-. 



If 7 = 0, then x = x(a) solves the equation x' = y/vx. Thus obviously x = c(a — a ) 2 and 
therefore 

/3 = </ -1 (a) = / = /3 0) d = -1/c. 

J c(a — «oJ a — «o 

The corresponding equation (|T3l) follows from (a — ao)(/3 — /?o) = ^- It is of the form 

(ui tl + u — a )(ux + ui — f3 ) = d. (24) 

Corollary of Theorem 2. If for a nonlinear chain ( Tl3p one of the conditions i) — in) is 
satisfied then the chain is of one of the form: 

1) U14 + u = c{u\ + ui + ci); 

aiMi^MiMiM + a 2 Ui^u + a3MiMi + a 4 = 0; 

3) (mi,i + u — a )(ui + Ui - Pq) = d. 

Let us study in details the list obtained. It is remarkable that it contains a new integrable 
example of equation ([!]). 

Theorem 3. The chain 2) in the Corollary of Theorem 2 admits a symmetry of the form 

u t = g(ui,u-i,u,ui,u-i) (25) 

if and only if the condition holds k := ^ = ±1. Under this condition the symmetry is of the 
form 

ku\ + u-\ _kux+U-i 

u t = Am- h [iu— = — . 

ku\ — m_i ku\ — U-i 

Chain 3) admits a symmetry of the form $25]) if and only if ocq = (3q and then it is reduced to 
the discrete potential KdV equation. 

Theorem can be proved by using technique developed in [5]. Actually it claims that the 
chain 2) with a 3 = ±a 2 passes the symmetry test. For the case k — 1 (or a% = a 2 ) the chain is 
a particular case of the Viallet equation (see [10J), this means that by a point transformation 
it can be reduced to one of the equations of ABS list [4j. Higher symmetries for this case are 
found in [6J. For some special choice of the parameters chain 2) with k = — 1 is reduced to the 
discrete sin-Gordon equation 

3(ui t iUiUiU — 1) + Ui^u — uiUi = (26) 

found by R. Hernandez Heredero and C. Scimiterna. Recently in [5] it was proved that equation 
( |26l) admits higher symmetries. To the best of our knowledge general case of the equation 2) 
with k = —1 has never been identified before as integrable one. 
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When a,i = and 0,3 = —0,2 it is reduced to the well-known discrete Liouville equation 

gUl,l+U _ gMl+Ml _|_ ^ 

having nontrivial integrals, i.e. functions solving the equations DF = F and DF = F, with 
F = e u ~ Ul + e U2 ~ Ul and F = e""" 1 + e" 2 "" 1 . 

The second part of the theorem claims that case 3) admits a symmetry of small order 
only if it coincides up to a point transformation with the discrete pkdv equation. The ques- 
tion is still open whether it admits any symmetry of a more complicated form, say u t = 
g(u2, Ui,U-x, u, Mi, U-i, U-2)- There are some technical difficulties with applications of the sym- 
metry approach to discrete models (see discussion in [5]). 

Turn back to the set G consisting of X, Y and their all multiple commutators. Assign two 
integers: order and degree to each element in G. Define the order of an element Z G G as 
a number of its factors X and Y minus one. For instance ord[X, Y] = 1, ord[X, [X, Y]] = 2 
and so on. Define the degree deg(Z) of Z as the exponent k in the expansion of the operator 
obtained by applying the automorphism (JHJ): DZD^ 1 = x k Z + where the tail is a linear 
combination of the elements with order less than ord(Z). Denote through Gij a subset of G 
containing elements with the order % and the degree j. Let Gi = Uj Gi,j be the union of all Gij 
with one and the same i. Evidently the set G^-i (as well as Gj -j+i) contains the only element 
^1,1-1 = ad l x (Y) up to the factor —1 (correspondingly, the only element Zj = ad % Y (X) up 
to the factor —1). Here the operator ad is defined as follows adx(Y) = [X,Y]. 

Theorem 4. Suppose that Z^k-\ ( or Zk,-k+i) is in the basis of the linear space Vk 3 Gt for 
any integer k: 3 < k < N, but Zn,n-i (respectively Z^-n+i) is linearly expressed through the 
other operators in Vjsr then function x = x(ui+U-.i) solves an equation of the form x' = e(x — l) 
with the constant coefficient e. 

The proof of the theorem is done exactly in the same way as the proof of Theorem 3 from 
|16j . but presented in the article for the readers' convenience. 

Lemma 3. For any integer k > 3 we have 

DZk+x^D' 1 = x k Z k+1<k + c k x x k ~ l Z k ^ k -\ H , (27) 

where c k > (but c k > for k > 3) and the tail contains a linear combination of the operators 
with the order less than k. 

Prove the Lemma 3 by induction. From the list of equations (1151) one gets for Z = Z 3 2 and 
Ri = Zifi the following representation 

'3 

— 12 / '2 " \ "' r 11 X / a 

DZ 3 2 D = x Z 3 2 + [x — 2xx )Z\ — qY — xqX, q = xx — 2x x H (28) 

x 

showing that the statement is true for the case k = 3. Suppose now that DZ k ^\D~ x = 
x k ~ 1 Z kjk ^ 1 + c k - 1 x'x k ~ 2 Z k -i :k _ 2 H and evaluate DZk+i^D' 1 : 

DZ k+1)k D~ l = [xX,x k ~ 1 Z kik ..i + c k -ix'x k ~ 2 Z k _ x ,k-2 H ] = x k Z k+hk + c k x x k ~ l Z k ^ x H 

where c k = c k -\ + (k — 1) > c k ^\ > 0. The proof is completed. 
Proof of the Theorem 4. Suppose that 

Z N ,N-l = a " Z v + b V Z V H > ( 29 ) 

ord(Z„)=N ord(Z,,)=N-l 
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where Z v and range the basis of V^j and the tail contains a linear combination of the 
operators of less order. Apply the automorphism (|HJ) to both sides of (|29|) : 



x 



( a " Z " + b V Z H H ) _ 2 ^'CiV-l^iV-l,7V-2 H 



ord(Z v )=N ord{Z^)=N-l 



Y, D(a v ){x k "Z v + •••)+ £ £>(M(^ M ^ + ---)- 




Collect the coefficients before and get 



x N ' l a v = x kv D(a v ), k u ^ N - 1. 



(30) 



Due to Lemma 2 functions a y and depend on x and its shifts. Moreover it follows from 
(130|) that a y cannot depend on x, x±i, x±2, ... at all. Therefore the only possibility is a v = 0. 
Compare now the coefficients before Zn-i,n-2 an d find: 



where b is the coefficient of Z^ -i i n t ne expansion (12"9"1) . A simple analysis of equation (|31|) 
shows that b is constant. Thus (I3T]) is equivalent to the equation x' = e{x — 1) with e = 6/cjv-i- 

Corollary of Theorem 4. T7ie case if) o/ i/ie Classification scheme never realized. 

Proof. Suppose on contrary that such a case is realized. Then at least one of the vector 
fields Zk t k-i or Z k _ k+1 should be linearly expressed through other elements of V k , otherwise 
A(k) > 2. Therefore due to Theorem 4 we have x' = e(x — 1) which corresponds to the case ii) 
A(0) = 1, A(l) < 2 and dimT = 3. The contradiction shows that our assumption is not true. 
The proof is completed. 

Let us summarize the results of reasonings above in the following theorem. 

Theorem 5. Suppose that an equation of the form (T3\) passes the algebraic test then it is 
of one of the form given in Corollary of Theorem 2. 

4 Conclusions. 

A classification scheme introduced recently in [16] is applied to quad-graph Klein-Gordon type 
equation. The list of equations passed the test contains along with the well-known integrable 
models some new example which passes also the symmetry test. 
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